All static spherically symmetric perfect fluid solutions of Einstein's Equations 



Kayll Lake 

Department of Physics, Queen's University, Kingston, Ontario, Canada, K7L 3N6 

(Dated: February 4, 2008) 

An algorithm based on the choice of a single monotone function (subject to boundary conditions) is 
presented which generates all regular static spherically symmetric perfect fluid solutions of Einstein's 
equations. For physically relevant solutions the generating functions must be restricted by non- 
trivial integral-differential inequalities. Nonetheless, the algorithm is demonstrated here by the 
construction of an infinite number of previously unknown physically interesting exact solutions. 
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Exact solutions of Einstein's field equations provide 
a route to the physical understanding (and discovery) 
of relativistic phenomena, a convenient basis from which 
perturbation methods can proceed and a check on numer- 
ical approximations. Here we look at static spherically 
symmetric perfect fluid solutions. Unfortunately, even 
for this simple type, very few solutions are in fact known, 
and of these few pass even elementary tests of physical 
relevance In this paper, an algorithm based on the 
choice of a single monotone function (subject to bound- 
ary conditions) is presented which generates all regular 
static spherically symmetric perfect fluid solutions of Ein- 
stein's equations. We are interested only in physically rel- 
evant solutions here and so the algorithm must be supple- 
mented by physical considerations These additional 
conditions limit the generating functions allowed by way 
of non-trivial integral-differential inequalities. The de- 
tails of how to choose physically relevant generating func- 
tions (beyond trial and error) are, at present, not known. 
Nonetheless, the robustness of the algorithm is demon- 
strated here by the construction of an infinite number 
of previously unknown physically interesting exact solu- 
tions. 



To set the notation, consider a spherically symmetric 
spacetime Jv{ Q 



First consider M. in "curvature" coordinates. 



ds 



M 



(1) 



where dil^ is the metric of a unit sphere {d9'^ + 
sill? {9)d(tP') and R = R{x^,x'^) where the coordinates 
on the Lorentzian two space S are labelled as and 
. Consider a flow (a congruence of unit timelike vec- 
tors m") tangent to an open region of S and write n" 
as the normal to u" in the tangent space of E. Both 
and n" are uniquely determined. We suppose that 
^ is generated by a fluid subject to the condition 
C^u^np = where is the Einstein tensor (see Q). 
Let G = G^, Gl EE G^u°'ufi and G2 = G^n'^np. In the 
static case it follows that the flow is shear free and that 



G + Gl = 3G2 



(2) 



is a necessary and sufficient condition for JQ) to represent 
a perfect fluid |^. 



dr'^ 



2m(r) 



r^dfl^ 



(3) 



Writing out ^ ^ we obtain an expression involving <I>(r) 
and m(r) with derivatives to order two in (f>(r) and to 
order one in m{r). Viewing as a differential equa- 
tion in ^(r), given m[r), we obtain a Riccati equation 
in the first derivative of $(r). However, viewing (O as 
a differential equation in m(r), given $(r), we obtain a 
linear equation of first order y^. As a consequence, we 
have the following algorithm for constructing all possible 
spherically symmetric perfect fluid solutions of Einstein's 
equations: 

Given <I>(r) (sufficiently smooth and subject to bound- 
ary conditions explained below) 



m(r) 



/b(r)e/°'''^'^''dr + C 

r a{r)dr 



(4) 



where 



^ 2r^($"(r) +$'(r)^) -3r$'(r) -3 
""^"■^ - r(r$'(r) -fl) 



and 



_ r(r(^"(r) + $ (r)^)~$ (r)) 
r$'(r) 1 

with = ^ and C a constant. The generating function 
associated with any known solution is of course immedi- 
ately obvious following the algorithm. 

Interior boundary conditions on <I>(r) are set by the re- 
quirement that all invariants polynomial in the Riemann 
tensor are finite at the origin. In this case there are but 
three independent invariants and these are expressed 
here in terms of the physical variables; the energy density 



P = 



Gl m (r) 



> 



(7) 



2 



and the isotropic pressure 



G2 {r){r — 2m{r)) - m{r) 



> 0. 



(8) 



Note that the inequalities in {Tj) and © are to be viewed 
as imposed restrictions on $(r). At the centre of sym- 
metry (r = 0) the regularity of the Ricci invariants re- 
quires that p(0) and p(0) be finite. The regularity of the 
Weyl invariant requires that m{r) is C'^ at r = with 
to(0) = m(0)' = m(0)" = and m(0)"' = 87rp(0) 9]. In 
summary, for a static spherically symmetric perfect fluid, 
finite p(0) and p(0) guarantees the regularity of all Rie- 
mann invariants at the centre of symmetry. $(0) is a fi- 
nite constant (set by the scale of t) and it follows from (jHJ 
that $'(0) = and $"(0) = f^(3p(0) + p(0)) > 0. Since 
p > and continuous and since p(0) > and finite it fol- 
lows from 101 that r > 2m(r) ^W\. With r > 2m (r) for 
r > it also follows from JSJl for p{r) > that $' (r) ^ 
for r > 0. As a result, the source function <i>(r) must 
be a monotone increasing function with a regular min- 
imum at r = 0. Exterior boundary conditions on <I>(r) 
exist only for isolated spheres and these conditions are 
set by junction conditions 0. The necessary and suffi- 
cient condition that A4 have a regular boundary surface 
with a Schwarzschild vacuum exterior at r = i? > is 
given by p{r — R) = 0. Setting m(r = R) = M it follows 
thatci>'(r = i?) = ,^(^^. 

Each source function $(r) which is a monotone increas- 
ing function with a regular minima at r = necessarily 
gives, via Q), a static spherically symmetric perfect fiuid 
solution of Einstein's equations which is regular at r = 0. 
Exact solutions, in the present context, can be viewed as 
those for which can be evaluated without recourse 
to numerical methods. The number of source functions 
$(r) for which can be evaluated exactly is infinite. It 
should be noted, however, that the generation of an ex- 
act solution does not necessarily mean that the equation 
p{r = R) = can be solved exactly. The algorithm pre- 
sented here is now demonstrated by the construction of 
an infinite number of previously unknown but physically 
interesting exact solutions of Einstein's equations. 



Let 



*W = o^ln(l + — ) 



(9) 



where N is an integer > 1 and a is a constant > 0. 
The function ^ is monotone increasing with a regular 
minimum at r = 0. With the source function @, Q 
can be evaluated exactly for any N. Whereas Q gen- 
erates a "class" of solutions, the metric (in particular 
m(r)) looks quite distinct, and the physical properties 
are quite distinct, for each value of N. Previously, only 
for N — 1,...,5 were solutions known, having been ar- 
rived at by various methods, and one solution which is 
the first term in the Taylor expansion of ^ . (These 



solutions, with N = 1,...,5, in fact constitute half of 
all the previously known physically interesting solutions 
(of this type) in curvature coordinates.) For N > 5 the 
solutions are acceptable on physical grounds and even 
exhibit a monotonically decreasing subluminal adiabatic 
sound speed 0. 

It is, perhaps, worth noting here that the foregoing 
discussion in curvature coordinates can be transformed 
directly into Bondi radiation coordinates 

Now consider "isotropic coordinates" 
dsl, = e2^('')(dr2 + T^dn^) ~ e2(*('-)-^m)dt2. (iq) 

Unlike curvature coordinates, the isotropic form H10() 
does not offer an immediate invaria nt p hysical interpre- 
tation of the functions ^'(r) or B{r) 15]. However, as we 
now show, the coordinates offer a simplified algorithm for 
constructing perfect ffuid solutions. Writing out |(5J we 
now obtain an expression involving ^'(r) and B{r) with 
derivatives to order two in ^'(r) and to order one in B(r). 
Viewing Q as a differential equation in *(r), given B{t), 
we again obtain a Riccati equation in the first derivative 
of ^(r). However, viewing |5J) as a differential equation 
in B(t), given vE'(r), we obtain an equation solvable sim- 
ply by quadrature. As a consequence, we have the fol- 
lowing simplified algorithm for constructing all possible 
spherically symmetric perfect fluid solutions of Einstein's 
equations in isotropic coordinates: 

Given ^'(r) (sufficiently smooth and subject to bound- 
ary conditions explained below) 



where 



c(r) 



S(r) = l'(r) + / c(r)fir-|-C 



V2 



(11) 



y^(,jf'(r))2„^"(i.) + ^'(r)/r (12) 



with e = ±1, = and C a constant. Recently, Rah- 
man and Visser [l^l have also presented an algorithm 
for constructing spherically symmetric perfect fluid solu- 
tions in isotropic coordinates. The source function ^'(r) 
used here is related to the source function z(r) used by 
Rahman and Visser as follows: 



*(r) 



rz(r) 



1 — z(r)r^ 



rdr. 



(13) 



The two algorithms differ fundamentally in the sense that 
only one integration is used in the present procedure as 
opposed to two distinct integrations used in the Rahman- 
Visser procedure. The Rahman- Visser procedure was 
motivated by the requirement that metric be manifestly 
real ab initio. The reality of the integral (|ll|l is discussed 
below. 
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Interior boundary conditions on ^'(r) are set exactly 
as in the case of curvature coordinates. We now have the 
energy density and pressure in the form 

(14) 

and 

P = g - ^i-B'ir)nr)HB'ir)r-2^) > 0. 

(15) 

\E'(0) is a finite constant (set by the scale of t) and it 
follows from that *'(0) = and from that 
B (0) = 0. With p{t) > it follows that the source func- 
tion ^'(r) must be a monotone increasing function with 
a regular minimum at r = and ^E* (0) = 47re'^'^^*'-'p(0). 
Exterior boundary conditions on ^'(r) are set as in cur- 
vature coordinates. Regularity of p{0) requires B (0) = 
and with p(r) > it follows that B{i) must be a mono- 
tone decreasing function with a regular maximum at 
r and B" (0) = -47re2^("V(0)- The limits -f < 
B (i) < guarantee the positivity of the effective gravi- 
tational mass. To examine the reality of the metric con- 
sider the function F{r) = (*'(r))^ - vl/"(r) -t- -^i' (t)/t. 
Now F(0) = 0, F'{0) = and F"(0) > so F{r) has 
a local minimum at r = 0. Now suppose that F(r) = 
for r > 0. Then condition ||2J) requires B — so we 
have already passed through a region with p < before 
the reality of the metric breaks down (in agreement with 
known theorems [Tol|'). 

In parallel to the algorithm in curvature coordinates, 
each source function ^'(r) which is a smooth monotone 
increasing function with a regular minima at r = neces- 
sarily gives, via a static spherically symmetric per- 
fect fluid solution of Einstein's equations which is regular 
at r = 0. Exact solutions are again those for which 
can be evaluated without recourse to numerical meth- 
ods. Physical considerations must guide the choice of 
^'(r). In isotropic coordinates, ratios of invariants and 
differential invariants can be obtained directly from the 
source function \E'(r) via differentiation. You do not need 
B{i) and in particular you do not need to integrate. For 
example, the functions p(r)/p(r) and p (r)/p (r) follow 
directly without integration. Of course, neither p(r) nor 
p(r) follow without integration. In curvature coordinates 
you cannot get these ratios without integration, starting 
from the source function $(r). 

To demonstrate the algorithm in isotropic coordinates 



let 

*(r) = aln44 (16) 
5(r) 

where a is a constant > 0. Of course it is not difficult 
to find functions /(r) and g{i:) so that Hlt)|) is monotone 
increasing with a regular minimum at r = 0. Nor indeed 
is it difficult to find such functions for which i?(r) can be 
evaluated exactly. For example, let g{i) = {6 + ei'^)'^ and 
/(r) = (5'»-f 7r^ with 6, e, 7 and C constants such that 6 > 
and (5^^''7 > (e. This class of solutions includes a num- 
ber of known solutions including the Schwarzschild inte- 
rior solution and the Rahman- Visser general quadratic 
ansatz. Any solution in isotropic coordinates can be im- 
mediately recovered and generalized following the algo- 
rithm presented 0| . 

An algorithm based on the choice of a single mono- 
tone function (subject to boundary conditions) has been 
presented which generates all regular static spherically 
symmetric perfect fluid solutions of Einstein's equations. 
In all cases the choice of generating function must be 
guided by physical considerations. These additional con- 
ditions limit the generating functions allowed by way of 
non-trivial integral-differential inequalities. The details 
of how to choose physically relevant generating func- 
tions (beyond trial and error) are, at present, not known. 
Moreover, the resultant equation of state is a byproduct 
of the algorithm and can not be set a priori. Despite 
these reservations, the algorithm has been demonstrated 
by the construction of an infinite number of previously 
unknown physically interesting exact solutions It 
is a curious fact of history that over half a century ago 
Max Wyman pointed out that the algorithm pre- 
sented here was possible and yet, despite the voluminous 
literature on the subject 0, the algorithm appears not 
have been followed up. 
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